Missing Quasiparticles and the Chemical Potential Puzzle in the Doping Evolution of 

the Cuprate Superconductors 
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The evolution of Ca2-i:Nai:Cu02Cl2 from Mott insulator to superconductor was studied using 
angle-resolved photoemission spectroscopy. By measuring both the excitations near the Fermi energy 
as well as non-bonding states, we tracked the doping dependence of the electronic structure and the 
chemical potential with unprecedented precision. Our work reveals failures in the conventional 
quasiparticle theory, including the broad lineshapes of the insulator and the apparently paradoxical 
shift of the chemical potential within the Mott gap. To resolve this, we develop a model where 
the quasiparticle is vanishingly small at half filling and grows upon doping, allowing us to unify 
properties such as the dispersion and Fermi wavevector with the behavior of the chemical potential. 

PACS numbers: 74.20.Rp, 74.25. Jb, 74.72.-h, 79.60.-i 



A central intellectual issue in the field of high- 
temperature superconductivity is how an antiferromag- 
netic insulator evolves into a superconductor. In prin- 
ciple, the ideal tool to address this problem is angle- 
resolved photoemission spectroscopy (ARPES), which 
can directly extract the single-particle excitations. De- 
spite the interest in this doping induced crossover, there 
continues to be a lack of experimental consensus, per- 
haps the most prominent example being the controversy 
over the chemical potential, /x. Over the past fifteen 
years, there have been conflicting claims of jj, either being 
pinned in mid-gap or shifting to the valence / conduction 
band upon carrier doping Q, |^ |M ^ IM 1^ ■ The inabil- 
ity of photoemission spectroscopy to provide a logically 
consistent understanding of this fundamental thermody- 
namic quantity has been a dramatic shortcoming in the 
field. 

In this paper, we present a new procedure to quan- 
tify /i with unprecedented precision by ARPES, while 
allowing simultaneous high resolution measurements on 
the low energy states. These measurements allow us to 
make major conceptual advances in addressing the dop- 
ing evolution. We find that the long standing confusion 
over /i stems from the manner in which quasiparticle-like 
(QP) excitations in the doped samples emerge from the 
unusually broad features in the undoped insulator. Our 
work reveals inconsistencies in the conventional frame- 
work that considers the main peak in the insulator spec- 
trum to represent the QP pole. On the one hand, we 
find that /i changes in a manner consistent with an ap- 
proximate rigid band shift; on the other hand, this shift 
appears to occur within the apparent Mott gap of the 
parent insulator. We show that this ostensible paradox 
can be naturally explained if one uses a model based on 
Franck-Condon-like broadening (FOB) where the quasi- 
particle residue, Z , is vanishingly small near half filling. 



This also reconciles existing puzzles regarding the insula- 
tor and the lightly doped compounds, and naturally ties 
the behavior of to low energy features such as the Fermi 
wavevector, kp, and the quasiparticle velocity v-p. 

Ca2_2:Naa;Cu02Cl2 is an ideal system to address the 
doping evolution of the cuprates. The stoichiometric par- 
ent compound, Ca2Cu02Cl2, is chemically stable and, 
along with its close variants, has served as the prototype 
for the undoped Mott insulator P, |^. Moreover, the 
system possesses a simple structure, with only a single 
Cu02 layer devoid of known superlattice modulations, 
structural distortions, or surface states, unlike the Bi- 
based cuprates, La2_£i;Srj;Cu04, or YBa2Cu307_5. The 
X = 0.10 and 0.12 samples had Tc's of 13 and 22 K, 
respectively (Tc,opt = 28 K), while the x = 0.05 compo- 
sition was non-superconducting, and were grown using 
a high pressure fiux method 7]. ARPES measurements 
were performed at Beamline 5-4 of the Stanford Syn- 
chrotron Radiation Laboratory with typical energy and 
angular resolutions of 13 meV and 0.3°, respectively, us- 
ing photon energies of 21.2 and 25.5 eV. Measurements 
were performed at 15 K, except for a; = 0, which was 
measured at T > 180 K. 

Previous quantitative studies of ^ have relied on core 
level spectroscopy. However, the precision of this method 
is limited and its interpretation can be very complicated. 
Since measuring /i is of paramount importance, we intro- 
duce a new approach which we believe to be more accu- 
rate and direct, and can be performed in parallel with 
ARPES on the near-Ep states. Our method utilizes de- 
localizcd, non-bonding O 2p states in the valence band, 
in particular, O 2pz at (0,0) and O 2pt^ at (tt, tt) shown in 
Figure^.. These were identified in earlier works and have 
no overlap with the Cu Sdx^-y^ orbital or Zhang- Rice 
singlet 0, [i3 . Because these well-defined peaks are rel- 
atively close to Ep and measured at a single wavevector. 
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FIG. 1: (a) Valence band spectra for x = 0, 0.05, 0.10, and 
0.12 compositions at k = (0,0) and {n,n). O 2pz and O 2p,r 
states are marked by triangles and circles, respectively, (b) 
Shifts of the O 2pz and O peaks shown on an expanded 
scale, (c) Doping dependence of fj,. 



we can treat these states as delocalized bands. The over- 
all shape of the valence band remains unchanged with 
doping, suggesting a rigid band shift on a gross scale 
(~ 10 eV). In Figuresnj) and^, we show the shift of the 
O 2pz and O peaks on an expanded scale with statis- 
tics collected from > 5 samples for each concentration. 
All data are referenced to the a; = composition, and 
we describe the methodology for determining /iq later 
in the text. The shift from these marker states yields 
A*o.05 = —0.20 eV, Ho.io = —0.28 eV, and ^0.12 — —0.33 
eV, all relative to /io, with a typical uncertainty of ±0.025 
eV. At finite a;, ^ = 1.8 ± 0.5 eV / hole, comparable to 
band structure 1.3 eV / hole) [ij]. 

These measurements reveal a fundamental failure of 
the traditional framework where the main peak of the 
insulator represents a quasiparticle pole, which we call 
the "coherent quasiparticle scenario" (CQS). Here, all 
energies above the peak maximum should fall within the 
Mott gap. As shown in Figure ^ /i shifts by an amount 
compatible with predictions from band structure calcu- 
lations. However, this shift appears to occur within the 
apparent Mott gap of the parent insulator - a logical in- 
consistency as there are no available states within the gap 
to shift into, as illustrated at the left of FigureEt^. While 
impurity-like states may form within the gap, in this pic- 
ture /i should not drop so rapidly with doping. In the 
remainder of this paper, we combine the measurements 
of jj, with high resolution studies of the near-Ep states. 
Our results suggest a model which borrows heavily from 
the concept of FCB, shown at the right of Figure [5^. 

To understand the failings of the CQS, the obvious 
starting point is the parent antiferromagnetic insulator. 
Early studies of Ca2Cu02Cl2 and Sr2Cu02Cl2 yielded 
broad peaks which exhibited a dispersion consistent with 



calculations for the extended t — J model P, 0, IT^ . 
These peaks were effectively interpreted as quasiparticle 
poles, in the context of the CQS. However, one crucial 
point that remained unresolved was the extreme width 
of these excitations. We address this as a critical weak- 
ness of the CQS, and use this as a starting point for 
constructing a new model. Understanding this lineshape 
is crucial, since many celebrated features such as the d- 
wave gap, the superconducting peak, dispersion anoma- 
lies, and the pseudogap [3], necessarily emerge from this 
starting point. Data taken at the top of the lower Hub- 
bard band, k — (7r/2,7r/2), is shown in Figure |2l3. In 
the CQS, one would expect the peak width, F, to be 
extremely narrow at the top of the band due to phase 
space constraints, analogous to excitations at the Fermi 
energy in a metal. Instead, F is comparable to the entire 
bandwidth 2 J ^ 350 meV, completely inconsistent with 
such a picture. Moreover, the width cannot be due to 
disorder, as the undoped system is stoichiometric, and 
adding chemical dopants results in sharper structures, as 
will be shown. For comparison, we also present spectra 
from Sr2Ru04 (thin black) exhibiting a nearly resolution- 
limited peak. Given that well-defined QP excitations 
can be observed by ARPES, we must confront the ori- 
gin of the broad peaks in Ca2Cu02Cl2. Moreover, in the 
CQS, the peak in Figurel^b should be well described by a 
spectral function yl(k, w) = , which 

should be approximately Lorentzian with a width dom- 
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FIG. 2: (a) Illustrations of the Coherent Quasiparticle Sce- 
nario (CQS) and the Franck-Condon Broadening (FCB) 
model, (b) Ca2Cu02Cl2 at k = {tt/2, 7r/2) with fits to a spec- 
tral function (dashed blue) and gaussian (red). A and B de- 
note the peak maximum and the onset of spectral weight, re- 
spectively. Comparison with Sr2Ru04 is shown (thin black). 
Up per inset shows photoemission spectra from H2, after Ref. 
(1^. (c) Dispersion of A and B along (0, 0) — (vr, tt), along with 
experimental values for (hues). 
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inated by an impurity scattering term, Fimp flj. A fit 
of ^(k, Lj) to the experimental data is shown and agrees 
poorly; to achieve even this, Fimp was assumed to be un- 
physically large (~ 300 meV), given that the material is 
stoichiometric and free of chemical dopants. 

In light of this failure of the CQS, we believe that an 
analogy to one of the simplest quantum systems, the H2 
molecule, may be enlightening. The H2 — > photoe- 
mission spectrum, shown in the upper inset of FigurelSb, 
exhibits FCB. Only the '0-0' peak (filled black) represents 
the final state with no excited vibrations and com- 
prises only ^ 10% of the intensity, while transitions to 
excited states with n = 1, 2, 3, and 4 vibrational quanta 
possess higher intensities than 0-0. In the solid state, 0-0 
alone would represent the QP or the coherent part of the 
spectral function, ^cohj whereas the excited states com- 
prise the incoherent part, Ainc- This behavior is redolent 
of polarons, and such models have been suggest ed in sys- 
tems where strong couplings are present [ijj uM ■ The 
low energy tail is suppressed exponentially, inconsistent 
with power law falloffs from -4(k, uj) spectral functions, 
but well-described by a Gaussian FCB envelope. 

Another unresolved issue is a large energy scale sepa- 
rating the peak from the experimental positions of /i in- 
side the Mott gap. For an insulator, /x is not well defined, 
and is pinned by surface defects and impurities and will 
vary between samples. However, the limits of this distri- 
bution are well defined, with a lower bound set by the 
QP pole at the top of the valence band. For this study, 
we identify two features : the peak maximum (A), and 
the onset of intensity (B), where B is determined from 
the first statistically significant signal above background 
(3(t). In Figure 13;, we show the dispersion of A and B 
along (0,0)-(7r, tt). While A qualitatively tracks the dis- 
persion of the t — J model, B disperses only weakly and 
has a large separation of > 450 meV from A. We present 
the distribution of /icxp from a large number of samples 
in Figure 13; and B clearly sets a lower bound for the dis- 
tribution of /icxp- This behavior suggests FCB where the 
true QP (B) is hidden within the tail of spectral inten- 
sity, and A is simply incoherent weight associated with 
shake-off excitations. For a; = 0, we reference A and 
the valence band features such that B is aligned to 0, 
and this is consistent with the lowest measured value for 
/icxp- This demarcates the upper bound for both the QP 
at half filling and /i(x = 0+), not the actual position of 
/i for X = Q, and provides us with the value for /xg which 
we have used throughout the paper. This model is also 
consistent with the temperature dependence of the line- 
shape, where a similar multiple initial/final state model 
was proposed 17]. At this stage, we cannot distinguish 
which interactions arc causing this broadening. Although 
numerical simulations of the t — J model predict that Z 
remains finite {Z ^ 0.2) j^^, some analytical calcula- 
tions have predicted that interactions with the antifer- 
romagnetic background cause Z ^ [mi. Another 
possibility is the coupling to the lattice, and our data 
bears some resemblance to that of lattice polaronic sys- 
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FIG. 3: (a-d) EDC spectra from a; = 0, 0.05, 0.10, and 
0.12 from (0.27r, 0.27r) to (0.727r, 0.727r) with hmnp positions 
marked by open symbols and kp shown in bold. Data are 
plotted on a relative energy scale referenced to the shift in fj, 
shown in Figure^. 

tems (ID Peierls, mang anites) [ll|ll|l3. We note that 
a recent calculation incorporating lattice effects in the 
t — J model has closely reproduced the ARPES spectra, 
including a vanishing QP peak (B) and a broad hump 
(A) which recovers the original t — J dispersion . 

In Figure 13 we show the doping evolution of the near- 
Ep energy distribution curves (EDCs), from (0.27r, 0.27r) 
to (0.727r, 0.727r). All data are plotted on a fixed energy 
scale relative to /io using the values determined in Fig- 
ure With doping, feature A evolves smoothly into a 
broad, high energy hump with a backfolded dispersion 
that qualitatively reflects the parent insulator (symbols), 
while fi shifts from B into the lower Hubbard band. It is 
now clear that /i does not fall immediately to A upon hole 
doping as expected in the CQS. Spectral weight develops 
at II, and a well-defined peak becomes visible for the x — 
0.10 and 0.12 compositions, comprising a coherent, low- 
energy band. The dispersion of the hump is summarized 
in Figure^, and was determined by tracking local max- 
ima in the EDCs; where the hump becomes less distinct, 
we also use the second derivative of the EDCs. We note 
that in our model, A no longer represents any precise 
physical quantity. We also track the dispersion of the 
lowest energy excitations (-0.05 eV < a; < Ep) from a mo- 
mentum distribution curve (MDC) analysis (lines). The 
collective behavior of the data reveals that A is roughly 
fixed at high energies (-450 meV) with doping, justifying 
our FCB model which, in some sense, decouples A from 
fi. The dispersion of the low-energy states reveals a re- 
markable universal behavior across doping levels where 
both the velocities of the QP dispersion (vp) and Fermi 
wavevectors (kp) virtually collapse onto a single straight 
line with a band velocity (1.8 eV-A) corresponding closely 
to the recently discovered "universal nodal velocity" [23| . 
This result also ties the chemical potential to the QP dis- 
persion to naturally explain how kp evolves with doping, 
by simply sliding fi down the true QP band and calcu- 



4 



0.0 



• 


x=0 




■ 


a:=0.05 




♦ 


j: = 0.10 




▼ 


.v = 0.12 ,' 


J ^lo.1o 



> 



^ -0.5 

LU 





0.2 0.4 0.6 

k (0,0)-{jr,Ji) 



0.8 



-0.4 0.0 
Energy (eV) 



FIG. 4: (a) Summary of hump (symbols) and MDC disper- 
sions (lines) from Figure|3 (b) Doping dependence of spectra 
from kp along with a schematic of the proposed distribution 
of coherent (blue) and incoherent (pink) spectral weight. 



lating Akp ^ Afi/vp, as shown in Figure^. We note 
that a number of different theoretical proposals have pre- 
dicted the emergence of sharp QP-like excitations from 
broad features, including dynamical mean- field theories 
|23j or "gossamer" superconductivity [2^ . 

A schematic cartoon overlaid on experimental data is 
shown in Figure ^J), illustrating the doping dependence 
of /i and the proposed transfer between coherent (blue) 
and incoherent (pink) spectral weight. This also explains 
the lack of any well-defined peak at kp for x = 0.05, ex- 
pected in Fermi liquid models, which can be naturally 



explained by the vast incoherent weight overwhelming 
any small coherent peak. This also clarifies whether 
in-gap states or a rigid band shift was the correct de- 
scription of the doping evolution, and it is now evident 
that neither are adequate in the FCB context. A similar 
emergence of QP excitations at low dopings was also ob- 
served in La2-2;Sr2;Cu04 ^2^, although one distinction 
between Ca2-a;Naa;Cu02Cl2 and La2_xSrxCu04 is the 
relative separation between A and Ep 

In conclusion, we have developed a phenomenological 
model based on high precision measurements of fi and 
detailed studies of the near-Ep states, providing us with 
the first globally consistent understanding of the doping 
evolution of the cuprates. This picture can be summa- 
rized as follows : i) At half filling, Z is vanishingly small 
in a manner reminiscent of Franck-Condon broadening. 
The true QP is found in the tail of spectral intensity, ap- 
proximately 450 meV above the peak position, ii) The 
previous misidentification of the peak maximum as the 
QP pole was at the root of the long standing confusion 
over fi. iii) With doping, spectral weight is transferred to 
the low-energy QP-like peak, iv) The shift of the chem- 
ical potential and kp is dictated by the band velocity of 
this faint QP band. We believe that the above picture 
provides a foundation for the origin of the quasiparticles 
upon doping and should be used as a guide to develop 
microscopic theories for high-Tc superconductivity. 
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